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Abstract. The effective long-time dynamics of solitary wave solutions of the 
nonlinear Schrodinger equation in the presence of rough nonlinear perturba- 
tions is rigorously studied. It is shown that, if the initial state is close to a 
slowly travelling soliton of the unperturbed NLS equation (in norm), then, 
over a long time scale, the true solution of the initial value problem will be 
close to a soliton whose center of mass dynamics is approximately determined 
by an effective potential that corresponds to the restriction of the nonlinear 
perturbation to the soliton manifold. 

1. Introduction 

In this paper, we rigorously study the long time dynamics of solitary wave 
solutions of the nonlinear Schrodinger equation in the presence of spatially rough 
nonlinear perturbations. Physically, nonlinear perturbations of the NLS equation 
arise in inhomogeneous nonlinear optical media, and in Bose-Einstein experi- 
ments where the scattering length, which determines the nonlinear coefficient in 
the Gross-Pitaevskii equation, is modified using Feshbach resonances, such as in 
experiments involving atomic soliton lasers, see [1] and references therein. 

Consider the initial value problem, 

(1) ^9iV^ = -A^-|^Ap> + e/(x,^), 

tp{t = o) = ^e H\R^X), 

where t G M denotes time, x G denotes a point in configuration space, with 
N > I, dt = the time derivative, A = YliLi ^ the iV-dimensional Laplacian, 

s G (0, ;|:), e G [0, 1) and / corresponds to a local nonlinear perturbation, which 
is a mapping on (complex) Sobolev spaces 

/ : H\R^,C) H-\R^,C), 

that is given by 

(2) /(x,^) = A(x)|^|2>, 

where s G [0, ^) (s G [0, oo) if = 1, 2,), and A G L°°(M^). 



^ Note that for a perturbation given by ([2]), (HJ is globally well-posed in H^, see for example 
Theorem 6.1.1 in [2]; and also [3]. 
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When e = 0, ([T]) admits travelhng sohton solutions, which are solutions of the 
form 

u{x,t) =e^^+*5'^-(^-'^)r/^(x-a), 
where v G is the velocity of the center of the soliton, a = ao+vt G denotes 
the center's position, 7 = 70 + G M is a phase, /x G M"*", and rj^ satisfies the 
nonlinear eigenvalue problem @ 

{-A + fi)7]^-r]l = 0. 

Starting with an initial condition which is close to a slowly travelling soliton, 
we want to study the effective dynamics of the true solution of the initial value 
problem if a rough nonlinear perturbation is switched on. We assume that 3cro = 
(ao, vo, 70, /io) G X X [0, 27i) x R+ and fixed a G (0, 1] such that 

(3) y - e^(^«+^^°(^-'^«»r/^„(x - ao)||^,i < c e^(i+"), 

N 
i=l 

for some finite constant c independent of e and a. We have the following theorem, 
which is the main result of this paper. 

Theorem 1. Consider the initial value problem (J\), and suppose that the nonlin- 
ear perturbation f is given by ([^, and the initial condition ip satisfies Then 
there exists eo > 0, such that for all 

eG[0,eo), /3g(0, a), i/ G (0, min(/3, a - 

there exists an absolute positive constant C , independent ofe, /3, and u, yet depen- 
dent on the initial condition, such that, for times <t < Ciy\loge\/e^^^^^~'^'^~°'\ 
the solution of the initial value problem (J\) is of the form 

^(x, t) = e*(^-("-")+^)(r/^(x - a) + w{x - a, t)), 

with 

where Oi = a — (3 — u>Q, and the parameters a, f , 7 and fj, satisfy the differential 
equations 

dta = v + 0{e^+^), 
dtv = -2^ V,K//,M(a) + 0(e^+"), 
1 . 2-Ns + 2s, 

-V 

4 2-Ns 



dfj = n + -v^ ^7^^ — 



9t/i = 0(e^+^), 



^See for example [5], Chapter 8, and references therein, for an overview. 
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where 
and 

We note that Vefj^^i corresponds to the restriction of the nonhnear perturba- 
tion to the sohton manifold, see Sect. O In other words, for initial conditions 
close enough to a slowly travelling soliton, and for small enough nonlinear per- 
turbations, the center of mass motion of the soliton is determined by Hamilton's 
(or Newton's) equations of motion for a point particle in an external potential 
corresponding to the restriction of the perturbation to the soliton manifold, up 
to small errors due to radiation damping and the extended nature of the soliton. 
The proof of this result relies on three basic ingredients: First, using a skew- 
orthogonal decomposition, Subsect. 12. 4[ the true solution of the NLS equation is 
decomposed into a part that belongs to the soliton manifold plus a fluctuation. 
By exploiting the group strucuture of the soliton manifold, Subsect. 12.51 and us- 
ing a Lyapunov-Schmidt mapping onto the tangent space of the soliton manifold, 
(155]) in Subsect. 13.1] it is shown that the dynamics of the part belonging to the 
soliton manifold is approximately determined by the restriction of the NLS equa- 
tion to the soliton manifold. Corollary [1] in Subsect. 13. 1[ As for the fluctuation, 
its norm is controlled over a long time scale using an approximate Lyapunov 
exponent. Proposition [2] in Subsect. 13. 2[ 

We now mention earlier results on solitary wave dynamics that are relevant to 
our analysis. In the last few years, there has been substantial progress in un- 
derstanding solitary wave dynamics of the NLS equation with time-independent 
potentials in the semi-classical limit, [U El El El El IS]- The latter analysis has 
been extended in [TD] to studying the interaction of a slow soliton of the cubic 
nonlinear Schrodinger equation with a delta impurity in one dimension, and in 
[TT] to studying the effective dynamics of solitons in time-dependent potentials 
in the space-adiabatic limit. The basic picture is that if the external potential 
changes spatially slowly compared to the size of the soliton (or when it is small 
in the case of a delta impurity), and if the initial condition is close to a solitary 
wave solution, then, over a long time scale (related to the external potential and 
the initial condition), the true solution of the initial value problem is close to a 
soliton whose center of mass dynamics is approximately that of a point particle 
moving in an effective external potential. The novelty in this paper is studying 
the effective dynamics of solitons in the presence of spatially rough and nonlinear 
perturbations, and also introducing in the analysis of soliton dynamics a useful 
Lyapunov-Schmidt mapping that is based on the group structure of the soliton 
manifold. We rely on extending different parts of the analysis in [8]- [11] to the 
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case at hand, and also on other important developments in the theory of nonlin- 
ear Schrodinger equations, [121 [131 [IH [U] ; see [21 E] for an overview. H We hope 
that along the way, we clarify basic general concepts that are helpful in the study 
of long-time dynamics of solitons. □ 

The organization of this paper is as follows. We first recall useful general 
properties of the nonlinear Schrodinger equation and the soliton manifold in Sect. 
[21 In Sect. [3l we prove the main result. Theorem [H 
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2. Some properties of the NLS equation and the soliton manifold 

In this section, we recall some properties of the nonlinear Schrodinger equation 
and the soliton manifold. We will use these properties in the following sections. 
In what follows, we denote by 

(4) g{u) := \u\^'u, ueH\R^,C), 

where s G (0, j^) appears in ([T]), and we let the functional G, F : H^(R^ , C) ^ M 
be such that G' = g and F' = /, where the prime corresponds to the Frechet 
derivative. The nonlinear Schrodinger equation is 

(5) zdttP = -Ai>-g{ij) + ef{x,tlj). 

2.1. Symplectic, Hamiltonian and Variational structure. The space H^(R^ , C) 
i7^(M-^,M^) as a real space, and it has a real inner product (Riemannian metric) 

(6) {u,v) := Re / dx uv, 



for u,v & H^lM.'^ ,C), where ■ stands for complex conjugation. It is equipped 
with a symplectic form 

(7) u{u,v) := Im / dx uv = {u,iv). 



■^We note that the problem of soliton dynamics in a slowly varying time- independent potential 
that was studied in [5] is revisited in [TB] for the special case of a cubic NLS equation in one 
dimension. 

'^We note that the analysis presented below can be directly extended to the case of more 
general nonlinearities, such a sum of local and Hartree nonlinearities, and to more general 
nonlinear perturbations, such as those with A slowly varying with respect to the size of the 
soliton, [17]. 

^The tangent space TH^ = H^. 



The Hamiltonian functional corresponding to the nonhnear Schrodinger equation 
© is 

(8) H^iij) ■.= 11 iVtPl^dx - G{ij) + eF{ij). 



2 _ 

Using the correspondence 

i7^(M^, C) < — ^ iJ^(M^, M) © if^(M^, M) 

■?/; < i> (Re?/^, Im-?/') 

< ^ J, 

where J := ( '^-i ] is the complex structure on if^(R'^,M^), the nonlinear 



-I 0^ 

Schrodinger equation ([5]) can be written as 

In the following, we denote H\R^) as either H\R^ , C) or H\R^ , M^). We note 
that since the perturbation is time independent, the Hamiltonian functional i^e 
defined in ([8]) is autonomous, and energy is conserved, i.e., 

(9) H,{^{t)) = H,{^{t = 0) = ip), teR. 
Moreover, is invariant under global gauge transformations, 

and the associated conserved Noether charge is the "mass" |^ 

(10) N{m) ■■=\j m)? = ^Mh, t e R. 

As mentioned in Sect. [T], when e = 0, the nonlinear Schrodinger equation ([S]) 
admits travelling soliton solutions that are of the form 

u{x,t) = e^^+^^''-(^-")r/^(x-a), 

where w e is the velocity of the center of the soliton, a = ao + vt E R 
denotes the center's position, 7 = 70 + G M is a phase, /i G M"*", and 77^ G 
L^(M^) n C^(M^) is positive, spherically symmetric and satisfies the nonlinear 
eigenvalue problem 

(11) (- A + /i)r/^ - r/J = 0, 
such that 

(12) |||a;|VllL2 + |||xnVr/^||U2 + \\\x\'^d^rijL^ < 00, 



^See for example [2], Chapter 6, for a proof of ^ and ([T 
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V/i G M^, and r/^ decays at oo as e~v^"^'". Furthermore, the "mass" of the sohton 



IS 



1 /" 2 l^E. 

(13) m(/i) = - dx 7]^ = 2 . 
Note that 

(14) a,m(/i) = (1 - y)^^"^'^ > 

since s G (0, j^), which imphes orbital stabihty, [HI [15]. We define 

(15) C,:=-A+f,-g'{r],), 

which is the Frechet derivative of the map (—A + — g{ip) evaluated at 
ry^. For all fi G M"*", the null space 

J^iC^,) = span{ir]f„ S^^.r/^, j = 1, ■ ■ ■ , A^}, 

see for example the Appendix in [8]. 

Orbital stability ((Hj) implies that is a local minimizer of ife=o(^) restricted 
to the balls Bm{n) '■= {"ip G : N{ip) = m(/i)}, where m(/i) is given in f[T3l) . see 
[T5] . They are critical points of the functional 



(16) S.iij) -.= 1 I dx (I V^p + M') - GW. 



2.2. Symmetries. When e = 0, the nonlinear Schrodinger equation ([5]) is in- 
variant under spatial translations 

Tf :^(x,t) ^^(x-a,t), a G R^, 

spatial rotations 

: ?/^(a;, t) ^ ^(i?-^a;, t), i? G 50(iV), 

time translations 

T^:i:{x,t) ^i:{x,t-r), 

gauge transformations 

T^:^{x,t)^e'^^{x,t), 7e[0,27r), 

and Galilean transformations (boosts) 

: tpix, t) e^'Xa;, t), ve M^. 

The conserved Noether quantities are the field momentum, angular momentum, 
energy, mass, and center of mass motion, 

^(-zV)V^, hixA-tV)^, ^ /"|V^A|2-G'(^), ^fm', /"^(a: + 2ztV)V. 



''See for example [5], Chapter 8, and references therein. 



Furthermore, when e = 0, ([5]) is invariant under complex conjugation 
and rescahng 

2.3. Soliton Manifold. We introduce the combined transformation 

T^, a = (a,i;,7,/i) G X X [0,27r) X M+, 
which is given by 

(17) T.^ := rjT^T^T^^iJ = e''~'^'^<^-^^+^^ fi^iVJ^ix - a), fit), 
where a, f G M^, 7 G [0, 27t) and /i G M"^. We define the soliton manifold as 

(18) Ms ■■= {Va = T^Vi -cr = {a, v, 7, //) G x x [0, 27r) x R+}, 

where rji satisfies (ITTil with = 1. The tangent space to the soliton manifold Ad^ 
at ?7o-o G A^s is given by 

'^vai.-^s = span{zt, Zg, Zb, Zs}, 

where 

^7-^7 '7o-o [7=0 ^''7(To 
2V^T^?7o-(, 1^=0 = ixrjf^g 

1 1 

h 

'2s 2 

In what follows, we denote 

(19) e, = -d^^, Ci+N = ixi, i = !,■■■ ,N, 



zt 



^9 



The soliton manifold Ais inherits a symplectic structure from [H^^uj). For 
cr = (a, t;, 7, /i) G X X [0, 27r) x M+, 

where is the orthogonal projection onto T^^TW^. It can be shown that VL~^ 
is invertible since d^m{fi) > 0, where the mass m(/i) = ^ f dx r]^ = , see 

[g. Exphcitly, 
(20) 



2N+2 



( ° 


— m(/i)l 





\ 


m(/i)l 




















m'(/i) 


V 





— m'(/i) 


/ 
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where e^?]^, k = 1, - ■ ■ , 2N + 2, are basis vectors of T^^Aig, 1 is the N xN identity 

matrix, and m(yu) = m'(/i) = 9^m(/i). Note that is related to ilj^^ 

by a similarity transformation. 

We note that the soliton rj^ breaks the translation and gauge symmetries of 
the nonlinear Schrodinger equation, which leads to associated zero modes of the 
Hessian C^, which is defined in f lT51) . Differentiating £'^{T^T^ri^) = with respect 
to a and 7 and setting the latter two to zero gives 

(21) C^zt = , C^zg = , 
while a direct computation gives 

(22) C^Zb = 2izt , Cf,Zs = izg. 



2.4. Skew-Orthogonal Decomposition. Consider the manifold A^'^ = {77^, a e 
So}, So = X X [0, 27r) x Jo, where Iq C I\dl and / C M+ is bounded. 
We define the 6 neighbourhood of in as 

Us:={tPeH\ inf \\^|J - < 6} . 

Then, for 6 small enough and for all ip G Us, there exists a unique (t{iIj) G 
Ci(f/5,S) such that 

^(V^ - V<T{i>), z) = {ip- ?7„(^,), J"^2) = 0, 

for all z e %f^^^^Ais- For a proof of this statement, we refer the reader to [S], see 
also [ig. 



2.5. Group Structure. The action of the combined transformation Tfj=(^a,v,-y,tJ.) 
defined in (fT7|) on elements of the soliton manifold Aig has a group structure 
Q /"Z, where Q corresponds to the semidirect product H^^+^ ix R"*", and H^^"*"^ is 
the real Heisenberg group in 2N + 1 dimensions. The quotient with Z is taken 
since 7 is defined modulo 27r. This group structure has been first noted in [TUj for 
the case = 1, but its generalization to higher dimensions is straightforward. 
The action of Q is given by 

(a', v\ 7', yu') ■ (a, V, 7, /i) = (a", v", 7", /j,"), 

where 

a" = + a', v" = J^'v + v' 
7" = 7 + 7' H -j=v ■ a' , and fi" = fi/i' . 

2y/ fi' 
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The anti-self adjoint operators {ej}j=i^... ^2N+2 defined in f fT9l) form the generators 
of the corresponding Lie algebra g. They satisfy the commutation relations 

(23) [cj, cj+n] = -e2N+iSij, 2, j = 1, ■ ■ ■ , A^, 
[ei,e2N+2] = ^ei, i = l,---,N, 

and the rest of the commutators are zero. Note that it follows from (flSjl that 
Aig is the orbit of Q/Z, and given z G %^^Ais, 3! G g such that z = Y^ri^. By 
exponentiating, the group element 

(24) r^={a,.,7,M) = e-'^-^-e''^e'^e'°s('^)(i+5-^-). 

3. Proof of the main result 

The proof of the main result is based on three ingredients, as discussed in the 
introduction. 

3.1. Reparametrized equations of motion. In this subsection, we use skew- 
orthogonal projection to decompose the solution of ([H) into a component that be- 
longs to the soliton manifold plus a fluctuation. We then use the group structure 
and a Lyapunov-Schmidt mapping to obtain effective equations for the parame- 
ters characterizing the component belonging to the soliton manifold. 

Given the initial value problem ([1]), let T be the maximum time such that the 
skew-orthogonal decomposition, defined in Subsect. 12.41 holds for tpit), <t <T, 
for some S > 0, i.e., there exists a unique a(t) such that 

ipit) = T^(j)?7i + w'{t), 

such that 

c^(w (t),FT^(t)?7i) = 0, 
for all y G g. In what follows, we denote r]i by rj. Let w{t) = T~^^w'{t). Then 

(25) z^(t) = T,(t)M(t), u{t)=7] + w{t). 

It follows from (123!) and the Baker- Campbell-Hausdorff formula, 

(26) e^Ye-^ = e'^'^^F, adx = [X, ■], 

for X, y G g, that T„YT^^ G g if F G g. Furthermore, it follows from translational 
invariance that uj{Ta-u,T„v) = fi~'~~uj{u,v), for m, t> G L^, and hence 

(27) uj{w{t),Yr^) = ^(T4i/;'(t),T;(J)(T.(,)rT-(,i))T.(,)7^) 

= ^^(^f)a;(u;'(t),FX{t)r/) = 0, 
VF G g, where Y' = T^it)YT,l G g. 
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We rewrite the nonlinear Schrodinger equation ([5]) as 

(28) dt{T^^t){v + wm = -iK{T,^t){r, + w{t)), 

where is defined in (jH]) and the prime stands for the Frechet derivative. 
Suppose that, for the interval [0,T], 

{a{t) = (a(t),^;(t),7(t),/i(t))}te[o,r] 

is a differentiable path in X X M X M+. Then, differentiating To-(t) with 
respect to t and using fl23l) and the Baker- Campbell-Hausdorff formula (|26|1 gives 

(29) dtT„(^t)=T,(t) X{a{t)) 
where X{a{t)) G g is given by 

N N . 

(30) X{a{t)) = ^ ciiy/JIei + ^ ^^e^+Tv + (7 - ^^)e2N+i + -e2N+2, 

i=l i=l V ^ ^ 

and the dot stands for dt. 
Using ([29]), it follows that 

dt{T^(t){v + ^W) = T.(t)[^(^^(t))r/ + X{a{t))w{t) + dtw{t)], 
and, together with (!28|) . this implies 

(31) 9i«;(t) = -X(a(t))r/ - X(a(t))«;(t) - zT-(|)i7:(T.(,)(r/ + ^(t)). 
Using ([H]) and ffT7|) in Sect. [21 a direct computation gives 



T-(,^)i/:(T.(,)(7^ + «;(t)) =^(r/ + w{t)) - /x[A(r/ + w{t)) + \r] + w{t)\''{r] + wit))] 



2 



- t^v ■ d^{7] + w{t)) + e/i»A(^^)|r/ + if(t)r'(r/ + w{t)). 

V A'" 

Substituting back in ([3T!) and using (fTTj) gives 
(32) 

where 
(33) 

N 9 



X{t) := -X(cr(t)) + ^ fiV/Iei + (;U - — )e27v+i 



N N . . 9 . 

E, . X /— —Vi , . a ■ V V fi 
[Vi - ai)y/fie, + ——ei+N + (-7 + ^3 -rJ^^'iN+i e2N+2, 

i=l i=l ^v/^ ^4 /X 

and 

£ = _A + 1 _ g'^r]) 
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the Hessian defined in f fTSj) . evaluated at 77, 

(34) g{r], w) = g{r] + w{t)) - g{ri) - g'{r])w{t), 

(35) f{x,r],w) = /i»A(x)[|?7 + tf;p^(?7 + w) - r/^^+^j. 
Note that it follows from (El), (jl), ([31) and ([35]) that 

(36) \\g{ri,w)\\H-i < C\\w\\jji 
and 

(37) \\f{x,v,w)\\H-^<C\\w\\Hi, 

uniformly in t G M, where C is a finite positive constant independent of e. 
We now introduce the Lyapunov- Schmidt mapping P : ^ g given by 

2N+2 

(38) P(-) := P„ 



a=l 



where e^, a = 1, ■ ■ ■ , 2N + 2, are given in (fT9l) in Sec. [21 and 

(39) Pi{-) := ^^^^^^j^^^YjUj{-,ei+NV) = -^{-^fii+NV), i = '^,---,N, 

Pi+N{-) ■■= —-^-—u;{-, CiT]) = u}{-, CiT]), i = l,--- ,N 
m(yU = 1) 

1 2s 
P2N^d-) := ;;^7(;^^(-'^2^+2r^) = ^3^-(-'^2^+2r;)' 

This mapping is well-defined since it follows from (fT2|) in Sec. [2l that X?7 G 
for all X = Yl^J^i"^ -^a^a £ g such that ||X|| := sup^^^i... 2N+2} \-^a\ < 00. Note 
that for X G g, 

2N+2 2N+2 

(40) P(Xr/) = Xo,Pp{eaV)ep = ^ X„5„;3e/3 = X. 

a,f3=l o,/3=l 

For t G [0,T], we introduce Y(t) G g that is given by 

(41) Y{t) := X{t) - e^^^P{^\C-^)r^''^^\ 

where X is given in f[33l) . The motivation for this choice of Y will become clear 
below. 

Proposition 1. Suppose satisfies Iji5\) such that ipit) G f/5, defined in Subsect. 
\2.4i for t G [0, T] and some 6 > 0, and let w, rj and Y be as defined above. Then, 
for \\w{t)\\H^ < c, c ^ (0, 1) independent of e, we have 

\\Ym<Cie\\wmm + Mt)rm), 
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where \\Y\\ := sup^g^j^ ... 2Ar+2} l^al ^'^^ C is a positive constant independent of 
t G [0, T] and e. 

Proof. It follows from ^ and ([39]) that 
(42) P{w{t)) = 0, 

and, since P is independent of t, P{dtw{t)) = dtP{w{t)) = 0. Together with fl32l) 
and ( l40l) . this implies 

(43) 

= P{dtw{t)) 

= Xit) + PiXit)wit)) - iiP{iCw{t)) - eiiip{i\(^^)Tf'^^) + iiP{icj{r], w)) 

- eP{if(^—^,r],w)). 

Since C is selfadjoint, it follows from ( l2Ti) and ( l22il in Sec. [2] and ( 127|) that 
uj{iCw{t), Caf]) = —uj{w{t), iCeaTj) = 0, a = 1, ■ ■ ■ , 2A^ + 2, 

and hence 

(44) P{iCw{t)) = 0. 
Now, dUD-dm) imply 

(45) Y{t) = -P{X{t)w{t)) - fiPmv. w)) + eP(^/(^, r/, ^)). 

We want to find an upper bound for the absolute value of the right hand side of 
(HSjl . In what follows, we denote by C a positive constant that is independent of 
e, and this constant changes from one line to another. 

It follows from Holder's inequality and (fT2l) in Sec. [2] that 

\uj{Xw,eo,r])\ < C\\X\\\\w\\l2, a = 1, ■ ■ ■ ,2N + 2, 

and hence 

(46) \\P{X{t)wm\<C\\X\\\\w{t)h.. 
It also follows from ([T2]) and that 

\uj{i~g{7],w{t)),e^7])\ < C\\w{t)\\l^, a = l,--- ,2iV + 2, 
and from ([12]) and (137]) that 



~ x + a^ r/, w{t)), e,r/)| < C\\wit) \\hi, a = 1, ■ ■ ■ , 2iV + 2. 



Therefore, we have 

m < rhn(fM\1... \\p(if(- 



(47) \\P{i~g{^Mm < C\\wmm, 11^ (^/(^, ^, < 
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Now, (!15])-(!17D give 

(48) < C{\\Ym\\w{t)h. + e\\w{t)\\H-. + 

and hence the claim of the proposition. □ 

We have the following corollary. 
Corollary 1. Suppose that the conditions of Proposition^ hold. Then, fort G 

dta = V + 0{€\\w{t)\\Hi + Mt)\\lr), 

dtv = -2^ VaVeff,^,{a) + 0{e\\w{t)\\H^ + H^i), 

^ dta ■ V 1 2 2 — A^s + 2s^^ , , s ^ , . 

+ O(6||./;(t)||Hi + imt)||^0, 
ai/i = 0(e||ii;(t)||Hi + lk(t)||?,i), 

where 
and 

Proof. It follows from fl41l) and Proposition [T] that 

(49) \\X{t) - e^,^P{zX{^W^')\\ < C{e\\w{t)\\H^ + \\w{t)\\U, 

for some positive constant C independent of e and t. 
Using ( fT2l) . it follows by integration by parts that 

(50) u;{tf{^,rj,w{t)),-d.,rj) = -J—l^'^d, [ dx\{^W'^\x), 



(51) 

,.~,x + a ,,,,,11 /I A^,i/",,,x + a, 25+2/ 



cj(z/(^— ,r/,u;(t)),(— + -x-ajr/) = ( )/i- / rfxA(^— )//'"+' (x) 



(52) dxA(-^)r/^^+^(x)x. 

Furthermore, since A and 77 are real, 

(53) uj(iX{^^)r]^'+\ 1X7]) = 

(54) ,^(,A(^±^)772m^,^) ^0. 
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The proof follows directly from ([33D, ([39]) and (|l9])-([5ll). □ 

3.2. Control of the fluctuation. In this subsection, we use an approximate 
Lyapunov functional and the coercivity property of the Hessian to obtain 
an explicit control on the norm of the fluctuation. Our approach is a slight 
modification of that used in [12], [13], [8]- [11]. Recall that it follows from the 
skew-orthogonal decomposition, Subsect. 12.41 that the solution of ([5]) can be 
rewritten as 

m = T,,,,iv,ix) + w'it, x)) ^ T'JT^T^irj^ix) + w'{x, t)), 

such that 

uj{w\Y7]^) = 0, VFGg. 

Here w' = T^w, where w appears in (125|1 in Subsect. 13.11 Let u' := rjfj, + w' . We 
define the Lyapunov functional 

(55) Ci^) := S^iu') + eF{^) - (S.iv,) + eFi^^^ix - a))), 

where £^ is defined in ( IT6!) . Sect. [21 and F' = /, where / is given in ([2[), Sect. [2l 
We proceed by estimating upper and lower bounds for C{ip). 

3.2.1. An upper bound for the Lyapunov functional. We have the following lemma. 

Lemma 1. Suppose ip satisfies such that ipit) G for t G [0,T] and some 
(5 > 0, and let u'.w'.rj^ as above. Then, for \\w'\\h^ < 1, there exists a constant 
C independent of e and t such that 

(56) \Cm))\ < C{yl + ey, + t{ezt + eyi + + vt)). 
where yt := sup^g^^t] \\w'{s)\\h^ and Zt := sup,g[o,j] \\v{s)\\. 
Proof. We have 



£fj.{u) ~ 2 J '^^l^'^'P + f^Wl'^ ~ G{u') 
= H,{u') + ]^^JL\\u'\\l,-eF{u') 



By translational symmetry. 
Furthermore, 



IT 



and hence 

(57) S,{u') = H.iij) + ^-i^-v' + /i)||^||i. - ■ V^) - eF(^). 



15 

Recall that since the perturbation / is time independent, energy is conserved, 

(58) dtH,{ij) = 0, 
while the rate of change of momentum is 

(59) dtiilj, -iV^) = e{f{x, ij), V^) + e{Vij, /(x, ^)). 

Formally, (1591) . which is a statement of a generalized Ehrenfest Theorem, follows 
from ([5]). To prove it, we introduce the regularizing operator Ih '■= (1 — hA)~^. 
For properties of 1^ as h 0, we refer the reader to Proposition 2.4.2 in [2]. 
Using ([5]), we have 

= hmiihi-Atfj - giij) + e/(x, ^)), hVtfj) 
/i— >o 

- (4V^, V4(-A^ - gi^) + e/(x, m} 
= e{(/(x,^),V^) + (VV^,/(V^))}. 

We also have 

where we have used in the first line the fact that 9t||^/'||i2 = (charge con- 
servation), and in the last equality the facts that ip = ei'"'^^~"'^^'^"'{r]fj_ + w'), 
{i'Vr]fj_,w') = {iw'yVr]^) = 0, and {ir]fj,,iw') = {iw',ir]^) = (which follow from 
skew-orthogonal decomposition). Furthermore, since rj^ is a minimizer of 

dtS^iv,.) = ^dtfi\\r]f,\\l2. 

It follows that 

dtCiij) = m^iij) - ^ev ■ {(/(^), V^) + (V^,/(^))} 



dMlvAh - ^9ta ■ VaF{ri^{x - a)) - edtixd^,F{'q^{x - a)) 



2 



1 dtv 



2 

+ -dtfi\\w'\\l2 - edta ■ VaF{ri^,{x - a)) - edtfxd^F{r]f,{x - a)) 
Together with f[T2]) . Corollary [T], and the fact that = ||T^_ii/;'||i/i < 



C||tf'||i:/i, it follows that 

(60) \dtCW\ < C{e\\v\\ + e\\w'\\l^ + e^u^'H^^i + 



\w 
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where C is a positive constant independent of e. Furthermore, by expanding 
^mo('7mo + ^o) around the minimizer ?7^g, we have 

(61) I^Mo(^w +^o) -^Mo(^Mo)l < C\\wXm- 
We also have 

(62) \F{ij)-F{r]^{x-a))\ = \F{r]^{x - a) + w'{t,x - a)) - F{r]^{x - a))\ 

< C\\w'\\hi. 

Now, fl60|) - (162|) together with the fundamental theorem of Calculus imply the 
claim of the lemma. □ 

3.2.2. A lower bound for the Lyapunov functional. In this subsection, we estimate 
a lower bound for C{ip). Let 

It follows from the coercivity property of that there exists a positive constant 

(63) p := inf {w, Cnw) > 0. 

We refer the reader to Appendix D in [8] for a proof of this statement. 

Lemma 2. Suppose ip satisfies ^ such that ipit) G f/^, t G [0,T], for some 
5 > 0, and let u',w',r]^ be as defined as above. Then there exists a positive 
constant C independent of e such that, for \\w'\\hi < 1, 

(64) \cim)\ > ^\\w'mi.-c{\\w'mi.+e\\w'{t)u.), 

where p appears in uniformly in t & [0,7"]- 

Proof. Expanding Sf^{u') around the minimizer 77^, we have 

(65) S,{r], + w') = SM + ^{w', C,w') + Rf{w'), 



where 



The claim of the lemma follows from 

\Rf\w')\<G\\wrH^, 

\F{4>) - F{r],{x - a))\ < C\\iv'\\h^. 
and the coercivity property ( l63i) . □ 
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3.2.3. Bound on the fluctuation. In this subsection, we combine both the upper 
and lower bounds on the Lyapunov functional and use a bootstrap argument to 
obtain a bound on the norm of the fluctuation. 

Proposition 2. Consider the initial value problem (|7]j, and suppose that the 
nonlinear perturbation f is given by and the initial condition ip satisfies 
Then, there exists eo > that depends on the initial condition, such that, for all 
e G [0,eo), P G (0,a) and v G (0,min(/3,a — /?)), there exists absolute positive 
constants C,C' , independent of e,f3, and v such that, for times 

\loge\ 



0<t<Cu 



max( 



(66) ||w/(t)||^i <CV+"-^-^ 

Proof. Suppose ip{t) G Us, t G [0,T], for some 5 > 0, and let u',w',rifj, be as 
defined as above. It follows from Lemmata [1] and [2] that, for yt < 1, there exists 
a constant C independent of e such that 

(67) y^ < C(yo' + eyt + yf + tiezt + ey"^ + e%t + yf)) 

where yt = sup^gjg.t] ll^'l-^)!!//! and Zt = sup^gjo.t] IK'l-^)!!- It ^^^o follows from 
Corollary [1] that 

(68) Zt<Zo + Ct{e + y^ + eyt). 

We choose /3 G (0, a), where a appears in ([3]), and we define 



(69) Tt := min( , ■ -^). 

e + yt + m ezt + eyf + e^yt + yt 

For t < n, dSTD and ([68]) imply that 

(70) yl<C'{yl + ey, + yl + e'^--P) 

<C"(^o + e"), 

for some positive constant C that is independent of e and (3. For e small enough, 
([3]) and the first inequality in ( 1701) imply that 

(71) yl<C{yl + e'+--P), 

for some positive constant C that is independent of e and (3. It follows that 

(72) y', < C{yl + 6^+"-^) 

<C(zo + e"), 

for an absolute positive constant C > 1 that is independent of e and /3, uniformly 
in t such that < t < r^. 

We now reiterate the above analysis. Consider the interval 

[0, T'] = [to, ti] U [ti, ta] U ■ ■ ■ U t„] C [0, T], 
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such that 

= to < ti < ■ ■ ■ < tn = T' , {U- ti-i) <Tt^, i = 1, - ■■ ,n. 

Let 

yi := sup \\w'{t)\\Hi, i = l,--- ,n. 
Zi := sup \\v{t)\\. 

t£[ti-i,ti] 

Note that from ([2]), 2/o ^ C"e5(^+") and zq < C"e", for some constant C indepen- 
dent of e. Iterating ( |72l) n times, we have 

n 

(73) yl < {J2 C^)C"e^+"-^ < c'^+^c'e^+^'-f' 

n 

We choose 

(74) e (0,min(/3,a-/5)). 
Given i/ and C, we choose n such that 

(75) C"+^ < €-". 
This imphes 

logC 

It follows from ([731) and ([75]) that 

(76) yl < C"el+"-^-^ 

(77) Zn < C'e"-^ 
We define 

(78) r :=^min( ^ ^ ' 



One can directly verify using ( l69l) . ( 1741) . ( 1761) and ( 178]) imply that that 

r' < Tf,, i = l, 

for e small enough. Therefore, 3eo > such that if e < eo, there exists absolute 
constants C and C independent of e, (3 and v such that 

uniformly for t E [0, uC\ loge|/e™°*^^~'^'^~°^], and hence the claim of the proposi- 
tion. □ 
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3.3. Proof of Theorem [T]. 

Proof. The proof foUows directly from CoroUary [T] and Proposition |21 □ 
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